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GROW-UP RATE AND REFINED ASYMPTOTICS FOR A TWO-DIMENSIONAL 
PATLAK-KELLER-SEGEL MODEL IN A DISK 

NIKOS I. KAVALLARIS AND PHILIPPE SOUPLET 



Abstract. We consider a special case of the Patlak-Keller-Segel system in a disc, which arises in the 

fl\ ' modelling of chemotaxis phenomena. For a critical value of the total mass, the solutions are known to be 

r^) ' global in time but with density becoming unbounded, leading to a phenomenon of mass-concentration 

, in infinite time. We establish the precise grow-up rate and obtain refined asymptotic estimates of 

'4) 1 the solutions. Unlike in most of the similar, recently studied, grow-up problems, the rate is neither 

polynomial nor exponential. In fact, the maximum of the density behaves like e^^' for large time. In 

particular, our study provides a rigorous proof of a behaviour suggested by Sire and Chavanis [Phys. 

Rev. E, 2002] on the basis of formal arguments. 
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1. Introduction 



■^ i 1.1. The complete Patlak-Keller-Segel model. Out of the many mathematical models that have 

^ ' been proposed to deal with particular aspects of chemotaxis, that proposed by Patlak in 1953 (cf. [iP] ) 

Q/ . and Keller and Segel in 1970 (cf. [21]) has received particular attention. The so called (two-dimensional) 

Patlak-Keller-Segel model consists of two equations, describing the evolution of the population density 

^^ . p(x, t) of bacteria, and the concentration c{x, t) of a chemical attracting substance, in a bounded domain 

' ' f2 C M^ and in a time interval [0, T] : 

'NT \ 

00 : % = y-{D,yp-xpyc), (1.1) 

0^ - Ac-ac + p. (1.2) 

More precisely the first equation describes the random (Brownian) diffusion of the population of cells, 
which is biased in the direction of a drift velocity, proportional to the gradient of the concentration of the 
chcmoattractant. The diffusion coefficient is denoted by Di > and the proportionality coefficient of the 
drift (mobility parameter) is denoted by x > 0. According to the second equation, the chcmoattractant, 
which is directly emitted by the cells, diffuses with a diffusion coefficient D2 = 1/0 > on the substrate, 
while is generated proportionally to the density of cells and at the same time is degraded with a rate 
equal to a/ 6 > 0. In order for system (|l.ip - (|1.2p to be well posed it should be supplemented with some 
initial conditions 

pix,0)^poix)>0, cix,0)^co{x)>0, (1.3) 

along with conditions on the boundary Qfi. A natural boundary condition, since it guarantees the con- 
servation of total mass, is the no-flux type condition for p, namely 

|^_p|£^0 on 917, (1.4) 

where v stands for the outer unit normal vector at dQ. As for c, a Dirichlct type boundary condition is 
assumed i.e. 

c = on dn, (1.5) 
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cf. [71 HH]. Note that the parabohc system (|l.ip - (|1.2p preserves the nonnegativity of the initial conditions, 
i.e. p, c > for t > 0, which is also expected to be true for the physical problem. For simplicity, Di,x 
arc considered to be constant and under suitable scaling can be taken Di = x = 1. 

In view of experimental facts, the coefficients and a are assumed to be small and a simplified form 
of the Patlak-Keller-Segel system is obtained (in fact this corresponds to the case when the diffusion and 
production of c are much faster than the dynamics of p and the degradation of c) . Namely, by considering 
the limiting case 9,a ^ 0+, the parabolic-parabolic system (|1.1[) - (|1.5[) is reduced to the elliptic-parabolic 
system 

^ = V-(Vp-pVc), xen,teiO,T), (1.6) 

-Ac = p, xeil, te (0,r), (1.7) 

-f-PjT = 0. xedn,te{0,T), (1.8) 

ov av 

c = 0, xe 917, f e (0,T), (1.9) 

p(a;,0) = /5o(a;)>0, x C Q.. (1.10) 

Note that in order for (|1.6p - (|1.10p to be well posed, only the initial data p(x, 0) = Pq{x) must be prescribed. 
Moreover, owing to the boundary condition (jl.Sp . the total (mass) population of cells is conserved, that is 

l|p(-,t)lli = l|Po||i=:A fori>0. 

For a more detailed analysis regarding the modelling as well as the behaviour of solutions of chemotaxis 
systems, see the review papers [MIEHIES] as well as the monograph [48]. Here, it should be noticed that 
system ()1.6|) - (|1.10p is also known as Smoluchowski-Poisson system and can describe the motion of the 
mean field of many self-gravitating particles [131 [H [51 [71 [IHl [321 [10] or that of polymer molecules [TB] . 
The behaviour of the solution to (|1.6p - (|1.10p . strongly depends on the parameter A. In fact, if A > Stt 
and Vl = S(0, i?), i? > 0, then solutions of (|1.6p - (|1.10p blow up in a finite time T*{po), that is 



lim \\p{-,t)\\Hi = \im \\p{-,t)\\LP = lim / p\ogp{ 

t^l t^l t^l Jj2 



a:, t)dx 



for every p > 1, see [9l Theorem 2(i)], [21 Theorem 2]. On the other hand for A < 87r, all solutions of 
system (|1.6p - p.l0p are global in time, cf. [7] Theorem 2 (iv)]. In the critical case A = Stt an infinite-time 
blow-up (grow-up) occurs, i.e. ||/3(-,i)||oo — » oo as i ^ oo, cf. [351 Theorem 3], [SI Proposition 3.2]. 
Finite or infinite-time blow-up can be accompanied by the occurrence of a 5— function formation in the 
blow-up set (this represents the trend of populations to concentrate to form sporae) and is known in 
the literature as chemotactic collapse. This phenomenon was conjectured by Childress and Percus |14| . 
Nanjundiah [38j . and was first verified, via matched asymptotics arguments, for a radially symmetric 
simplified Patlak-Keller-Segel system in [5S]. A result regarding the infinite-time Dirac mass formation 
for p can be found in [39] , where some characterization of grow- up (mass-concentration) points together 
with more grow-up results for different types of boundary conditions are also obtained, [39[ Theorem 
2 & Theorem 3]. For blow-up results concerning a variation of system ()1.6p - (|1.10p but with Neumann 
boundary conditions for both p and c, see |30l [35l [H [36l [44] . 

1.2. The simplified Patlak-Keller-Segel model in the radial case. In the case when H. is the ball 
B{0,R), R > 0, and the initial data po{x) = po{r) is radially symmetric, the solution of system (|1.6p - 
(jl.lOp is radially symmetric, i.e.: p{x,t) ~ p{r,t), with r = |a;|. In this case the elliptic-parabolic system 
(|1.6p - (|1.10p can be greatly simplified. Namely, by introducing the cumulative mass distribution 



Q{r,t) := / p{x,t)dx ~27r / sp{s,t)ds, 

J B(Q,r) Jo 

which is equal to the mass contained in the sphere B{0, r), the system reduces to a single equation 

Qt = Qrr - -Qr + 7;—QQr, < r < i?, t > 0, (1.11) 

r 2nr 

g(0,f) = 0, Q(i?,t)=A. (1.12) 
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By the definition of Q, the function 

Qir,0)^Qo{r) (1.13) 

is positive nondecreasing and satisfies the compatibihty conditions Qa{0) = and Qq{R) = A. 

As mentioned in [H [7] the formulation (|1.11[) - (|1.13[) aUows the consideration of some initial data for 
the density p which could be either unbounded or singular (such as measures), a case that seems rather 
realistic. This means that the initial data Qq for problem p.lip - p.l3p could have unbounded derivatives 
Qo.r or even be discontinuous. Moreover, using formulation p.lip - (jl.l3p we have the comparison principle 
at hand, which is not available for system (|1.6|) - (|1.10|1 . Due to the scaling properties of equation (jl.lip . 
we can assume without loss of generality that problem (|l.ll[) - p.l3[) is posed in the unit ball -8(0, 1). 
(Indeed, it is easily seen that if Q{r, t) is a solution of (|l.lip - (|1.13p then Q{Rr, R^t) is also a solution.) 

Using the new variable x = r^ (no confusion with the original variable x in p.6p - (|1.10p should occur) 
and defining N{x, t) = Q{r, t), we are led to the problem 

' t>0, (1.14) 



Nt = 


= AxN^^ + -NN^, < a; < 1 


N{0,t) = 


= 0, A^(l,t)-A, 


N{x,0) = 


= ^o(a;), < a; < 1. 



(1.15) 
(1.16) 

Note that (|1.14p differs from the Burgers equation only by the variable coefficient x in the diffusion term. 
The above problem, although it degenerates at a; = 0, may be handled more easily than (|l.lip - p.l3p 
since it contains less terms, and at the same time does not have any singular coefficients in the first order 
terms. 

As is expected, the behaviour of the solution of problem (|1.14l) - (|1.16p (which is well defined for suitable 
initial data) depends on A. For A > Stt the solution N ceases to exist in a finite time T*, more precisely 
the boundary condition A^(0, f) = is no longer fulfilled at t = T*. Moreover a "gradient blow-up" occurs 
at i = T*, in the sense that N^iO, t) —* oo as t s- T* and the density p becomes also unbounded at time 
T*, cf. [H Theorem 2(i)]. On the other hand, for < A < Stt and any (admissible) initial data there is a 
unique global-in-time solution N e C([0,oo);i^(0, 1)) n C^'^((0, 1) x (0,oo)). Furthermore, N converges 
to the unique steady state solution : 

N{-,t)^ Nd = 8TT—^, asi^oo (1.17) 

X + a 

in LP(ri), p> 1, and even in L°°{n) provided that supj>Q lA^^^loo < oo, where d = ^ — 1 > 0, cf. 3 . In 
this case the rate of the L^-convergence of N{-,t) to Nd is shown to be exponential. 

In the borderline case A = Stt, the situation is still different and the problem exhibits a typical critical 
behaviour. Namely, it is proved in [S] that there exists a global-in-time solution N, which converges to 
the "singular" steady state N{x) = Stt (note that N does not satisfy the boundary condition at a; = 0). 
Actually, as was proven in [391 Theorem 3], an infinite-time Dirac mass formation at the origin r = 
of the ball occurs in this case. However, neither an estimation of the grow-up rate nor the asymptotic 
profile of the grow- up are provided in |39| . On the other hand the authors in [SI Proposition 3.2] obtain 
the decay estimate 

||A^(-,t)-8^||Li <^ for t>l. (1.18) 

Estimate (|1.18p seems to be far from optimal since formal asymptotics performed in g2] suggest a 
temporal decay estimate of order 

||iV(-,<)-87r||i,i «0(e"^, as i -^ oo. (1.19) 

Remarks 1.1. (a) The (nonradial) Patlak-Keller-Segel system has also been studied in the whole 
plane K^. Again, the behaviour of solutions depends on the initial mass of the system and a dichotomy is 
found [HllTTj. More precisely, assuming < (1+ |a;p)/9o and pologPo & L^, there exists a critical value 
of the mass Nc := Sir such that if < ||/3o||i < A^c (subcritical case) only global- in-time solutions exist, 
while if IIpoIIi > A'c (supercritical case) the solutions blow up in finite time [TTJUT]. Moreover, in the 
subcritical case, solutions converge to a self-similar profile as t ^ oo [S1[TT]. Finally, for the critical case 
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N = Nc, which was studied in [10] . the solution is global- in-time and grows up as a Dirac mass at the 
centre of mass as f — > oo. 

(b) The only previous mathematical study of grow-up rates for a system of Patlak-Keller-Segel type 
concerns high dimensions, namely n > 11, and was performed recently in |45| . There, some radial global 
unbounded solutions were constructed in a ball and an infinite sequence of polynomial grow-up rates was 
obtained (for a suitable sequence of initial data). On the contrary, our results in the present paper for 
n = 2 exhibit a grow-up rate independent of the initial data, 

(c) Concerning the parabolic-parabolic Patlak-Keller-Segel system, in a bounded domain with Neu- 
mann conditions, interesting results can be found in [27] and [19], respectively on the asymptotics of 
finite-time blow-up, and on the convergence of bounded solutions. 

(d) Problem p.l4p - p.l6p is in fact independent of the boundary condition p.9p for c. Of course the 
boundary condition for c has to be taken into account if one wants to determine c, and not only p, from N. 

(e) In the case where the diffusion of cells is very slow as compared to the diffusion of the chemoattrac- 
tant, i.e. when Di ^ D2 in (|l.ip - (|1.2p . the complete system (|l.ip - (|1.5p is reduced to a single but non-local 
equation, cf. [51| . The global existence (subcritical case) and the finite-time blow-up (supercritical case) 
of solutions of the derived non-local equation in the two-dimensional case is studied in [31] . 

1.3. Heuristic description of the results and methods. Our aim is to prove rigorously the decay 
rate p.l9p . as well as to provide a refined asymptotic profile for N{x,t) as i — > 00. To normalize the 
constants arising in the calculations we shall work with the equivalent problem 

(1.20) 
(1.21) 
(1.22) 

which is obtained from p.l4p - (|1.16p by setting N{x,t) = 8TTu{x,4:t). 

It is clear that the stationary part xu^^ + 2u.Ui, = of the parabolic equation p.20p is invariant under 
the rescaling u{x) 1-^ u{kx) (fc > 0). Moreover, the steady state solution for ^ < 1 is given by 

Ua{x) = 1 — - = -— - 

ax + 1 aa; -|- 1 

for a = U'^{Q) — jl-? > 0. Rewriting (|1.17p in terms of u we have 

u(x,t) -^■ Ua{x), as t — > 00. 

Also, observe that Ua{x) converges in a monotone increasing way to the "singular" steady state U = I 
as a — > 00 (^ ^ 1). 

Motivated by the above considerations, we shall look for upper and lower solutions of problem (|1.20p - 
p.22p . which are perturbations of a moving family of steady states 

Uait){x), (1.23) 

the perturbation being defined in terms of the self-similar variable y ~ a{t)x. Here, a is a function of time, 
diverging to infinity, which is a priori unknown and will be eventually identified by a suitable "matching" 
procedure (see Remark 1.3). Note that such a form was used in [46| to construct "approximate solutions" 
leading to the formal asymptotics mentioned above. However, the expansions in [46] contained only a 
correction term at the first order, while those that we here construct involve first and second order terms 
(see Lemmas 14.11 and 14. 2p . This seems necessary to obtain (rigorous) upper and lower solutions living 
"close" to the actual solution, and eventually providing us with the desired decay rate, as well as with a 
good description of the asymptotic profile. 

As a consequence of this construction we shall derive an infinite-time boundary gradient grow- up result, 
with the grow-up rate 

M^(0,i) = A(<)(l + 0(t"^/^logi)), asi^oo, where ^(i) = exp - + V2i . (1.24) 



Ut = 


= xuxx + 'i.uux, < a; < 1, f > 0, 


u{Q,t) -- 


= 0, u{\,t)^i:^—, t>0. 


u{x,0) = 


= Ui)[x), < X < 1, 
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In terms of the solution of the original system (|1.6p - (ll.l0p . p.24p gives an estimate of the central density 
of bacteria, since p{0,t) = 8ux{0,4:t)- Note that (|1.24[) was also predicted by formal arguments in [46] , 
At the same time, we show the C^ regularity of u up to the boundary for all finite time intervals, hence 
ruling out the possibility of sup^-grg ii Ux{x,t) blowing-up in finite time. This problem was left open in [S]. 
Moreover, we obtain a precise asymptotic expansion of the solution (see formula (|2.8p in Theorem l2.ip . It 
expresses the solution as the sum of a quasi-stationary profile (cf. (|1.23p ) and of a correction term which 
becomes significant only for x bounded away from 0. As a consequence, we obtain the decay 



h(-,<)-l|Ui(o,i) = (l + 0{t-'/^logt))V2ie^v 



l-VTt 



as t — !• oo, 



again in accordance with the predictions in 



Remark 1.2. The study of unbounded global classical solutions of super linear parabolic problems and 
their asymptotic behaviour has recently attracted substantial mathematical interest. Particular effort has 
been devoted to the reaction-diffusion equation 

ut- Au^ uP, (1.25) 



where such solutions are known to exist for suitable p > 1. Let us mention the works [33l [181 HI] for the 
Dirichlet problem in a ball, and [42l [20l [2T] [34] for the Cauchy problem. See also [43l Sections 22 and 
29] for related questions. The case of the Frank-Kamenetskii equation (with nonlinearity e" instead of 
uP) is also studied in [TS]. As for the diffusive Hamilton- Jacobi equation 

Ut ~ Uxx = luxl'' , (1-26) 

with p > 2, results of this kind can be found in |47j . A common feature in all these examples is the 
stabilization of the solution to a singular steady-state, the growing-up quantity being ||w(t)||oo (resp., 
||Ma;(i)||oo) for equation (|1.25p (resp., (|1.26p ). It has to be noted that the grow-up rates, as i ^ oo, 
usually behave either like e^* or t'^. The only known exception (see [52]) seems to be the case of equation 
(|1.25[) with zero boundary conditions, in a ball of W^ with critical Sobolev exponent (p = 3). In this 
situation, unlike in spatial dimensions n ^ 4, there holds log ||u(t)||tx) ~ 2v^, thus leading to a similar 
rate as in our problem. 

Remark 1.3. Let us point out that our determination of the grow-up rate for problem (|1.20p - (|1.22[) is 
achieved through the matching of (sub- or super-) solutions with the imposed boundary condition at the 
right end-point a: = 1, an idea also present in [22]. This is different from what is done in [T8 [ [20 l [21 ] |47] . 
where the grow-up rate is determined by the matching between inner and outer (sub-/super-) solutions. In 
those works, the inner solution corresponds to a self-similar, quasi-stationary evolution along a continuum 
of regular steady states (similar to (|1.23p ). but the outer solution is obtained by a linearization around 
the singular steady state. Here, on the contrary, the behaviour of u in the inner and outer regions is 
unified through a single self-similar variable y = a{t)x (cf. formulae (|4.ip and (|4.27p below). Moreover 
we point out that, in our case, a linearization around the "singular" steady state U = 1 would not give 
the desired grow- up rate (jl.24p , but only a non-optimal exponential upper bound based on an associated 
eigenvalue problem. 

The paper is organized as follows. In Section 2 we state the main results. Section 3 contains a 
number of preliminary results: In Subsection 3.1, we recall the basic facts concerning local existence and 
comparison. In Subsections 3.2 and 3.3, we show that a control on the slope at a; = is enough to prevent 
gradient blow-up (Lemma l3.2p . and we obtain preliminary estimates of solutions for small time, which will 
be useful in Section 4 to initialize the comparison with the main sub-/supersolutions (Lemmas 13.31 and 
13. 4p . Subsection 3.4 is devoted to the study of a second order ordinary differential operator which plays a 
key role in the subsequent construction of sub-/supcrsolutions. The proofs of the main results are given 
in Section 4: Subsections 4.1 and 4.2 are devoted to the construction of the main sub- and supersolutions 
respectively; The proofs of Theorem 12.11 and Corollarv l2.2l arc finally completed in Subsection 4.3. 
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2. Main results 



Consider the problem 

ut — xuxx = "^uux, < a; < 1, t > 0, (2-1) 

u(0,i) = 0, t>0, (2.2) 

u{l,t) = 1, t>0, (2.3) 

u(x,0) = Uo{x), 0<a;<l. (2.4) 

Concerning the initial data, we assume that 

uq G C([0, 1]), uo(0) — 0, uo(l) = 1, "0 is nondecreasing, (2.5) 

and that 

uo{x) < Kx, < a; < 1, for some K >1. (2.6) 

Problem (P?T|) - (P^ admits a unique global solution, with u e C([0, 1] x [0,oo)), u G C^'^{{Q,1] x (0,oo)) 
(see [5] and Section 3.1 below). Moreover it was shown in [5] that < u < 1, m is nondecreasing in x, 
and 

lim u(x,t) ~ 1, uniformly for x in compact subsets of (0, 1]. (2-7) 

t — ^oo 

Our main results are the following: 



Theorem 2.1. Let uo satisfy 112. 5\) and 112. 6\) and denote by u the global solution of problem Ii2.1\) -t2. 
Then there holds 

l-x + 0(t-^/'^\ogt) 
1 — u{x,t) = —r— , uniformly in [0, 1], as t ^ oo, (2-8) 



with 

'5 



2 ■ ■'' 



(2.9) 



A{t) = exp 

Moreover, we have the regularity property 

w^eC([0,l] X (0,oo)) (2.10) 

and the estimate 

Ux{0,t) ^ A{t) (l + 0{t'^^'^ logt)Y ast^oo. (2.11) 

Corollary 2.2. Let uq satisfy \2.5]) and 112.6]) . denote by u the global solution of problem i2.1]) - (K4^ , 
and let A{t) be given by 112. 9\) . 

(i) The solution u satisfies the inner layer expansion (quasi- stationary behaviour) : 

1 — u{x^ t) = — - — , uniformly in any region x < o(l), as i — > cx). 

1 + A(t)x 



(a) We have the L^ -decay rate 



ii-,t) - l|Ui(o,i) = (l + 0(t-i/2 logf)) V2i exp 



-|-^/S 



as t ^ oo. 



Mx>0, 


< a; < 1, t > 0, 


Ut — XUxx ~ 


2uux, 0<a;<l, t > 


lim u(x, t) ~ 


0, for a.e. t > 0, 


u{l,t) = 


1, <>0, 


m(-,0) = 


uo inLi(0, 1). 
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3. Preliminaries 

3.1. Local existence and comparison principle. By [51 Theorem 2.1], we know that for any mq 
satisfying (|2.5p . problem (|2.ip - (|2.4p admits a global solution u in the following sense: 

ueC{[0,^);L\0,l)), (3.1) 

uG C2^i((0,l] X (0,oo)), (3.2) 

(3.3) 
and 

(3.4) 
(3.5) 

(3.6) 

(3.7) 

Note that the solution is obtained in [5] as limit of solutions of regularized problems, where xu^x is 
replaced by {x + e)uxx, e > 0. Also, for each T > 0, w is the unique local solution of p.ip - p.7p on (0, T). 
If, moreover, uq satisfies (|2.6p . then 

ueC([0,l] X [0,oo)), 

see also [H Theorem 1 (i)]. The continuity for t > and x ~Q follows from Propositions 2.4 and 2.5]. 
For t = Q and x G [0, 1], the continuity can be established by comparison with simple barrier functions. 

The following proposition provides a comparison principle suitable to our needs. 
Proposition 3.1. Let r > and the functions u,v satisfy the following regularity conditions: 

u,ue C([0,r);Li(0,l)), (3.8) 

u,weC2'^((0,l] X (0,r)), (3.9) 

Ux,VxeLl^{[0,l]x{0,r)), (3.10) 

u,t;eL,-([0,l]x(0,r)). (3.11) 

Assume that 

Ut — xuxx — 2uUx < vt — xvxx — "^vvx, 0<a:;<l, 0<t<T, (3.12) 

lim u(x,t) < lim v{x,t), for a.e. tG (Oi''")j (3.13) 

x-*0+ 2-->0+ 

u{l,t) < v{l,t), 0<t<T, (3.14) 

w(-,0) < v{-,0) a.e. in {0,1). (3.15) 

Then u < v in (0,1) x (0,t). 

Observe that (|3.10p implies that the limits in (|3.13p exist for a.e. t G (0,r). Note also that for 
any uq satisfying (|2.5p . the solution u of problem (|3.ip - p.7p satisfies conditions p.lOp and (|3.1ip . as a 
consequence of p.3p . (|3.5p and (|3.6p . 



Proof of Provosition \3. 11 It is a modification of the stability proof in [5} Theorem 3.1]. Let z = u — v. 
By (12321), we have 

Zt < xzxx + '2uUx ~ 2vVx ~ dx{xzx + z{u + V — 1)), 0<a;<l, < t < t. (3.16) 

For 5 G (0, 1) we define the following C^ (and piecewise C^) convex approximations of the function 
s —>■ s-^- = max(s, 0): 

0, if — oo < s < S, 
?i5(s) = <( (2(5)-i(s-(5)2, iid<s<2S, 
3(5/2, if 2(5 < s < oo. 
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Fix < ti < <2 < T and S G (0, 1). Then, for any e G (0, 1), multiplying p.l6p by (/>^(z), integrating 
by parts and using (|3.1ip . (|3.14[) . < (j)'g < 1 and (j)'g > 0, we obtain 

cj)s[z{x,t2))dx - I (f>s{z{x,ti))dx 



(l)g{z)dx{xzx + z{u + V — 1)) dxdt 



(j)'g{z)zt dxdt < 



t2 



ti 



4>'s{z)(xZj: + z{u + V - 1)) 



dt 



<j>'g{z)(xZr^ + z{u + f — l))^:!- dxdt 



4>g{z){u + f — l)zZx dxdt 



Je+K, 



< e f'\zxie,t)\dt + C f\'sizie,t))dt~ f' 
Jti Jti Jti 

where C depends on ti,t2 but is independent of e (and 5). Owing to (j3.10p . there exists a sequence 
e„ — > 0+ such that hm„^oo -fe„ = 0. Next, since hm„^oo ^^(^(sn, t)) = for a.e. t e (^1,^2) due to (jS.lSp 
and the definition of (f)s , we deduce that lim,, 



1 .1 

(j)s{z{x,t2))dx- (j)s{z{x,ti))dx < - 
Jo Jt 



n — >oo 



Je„ = by dominated convergence. Consequently, 

z4>'g{z)\\zx\dxdt. 



1 l-t2 /•! 



(f>^{z)(u+v~l)zZxdxdt < C 
Jt 



Now, observe that lim5_>o 0<5(s) = s+ and lim^^o ^(l^sis) = for each s £ M. Using < 4>s{s) < s+, 
()3.10p and (|3.1ip . we may pass to the limit S —^ hy dominated convergence in the preceding inequality 
and we obtain 



Zj^.[x,t2) dx — I z^{x,ti) dx <Q. 
/o Jo 

Letting ti -^ 0+ and using p.Sp and p.lSp . wc conclude that Jp zj^{x, t^) dx — for all ^2 G (0, r), hence 
u < u in (0, 1) X (0, r). D 

3.2. Sufficient condition for C^ regularity. As noted in, e.g., [12l Section 2.2] and in [23l Section 2], 
by means of the transformation 

Q 1 /" 

w{r, t) := ^u{r^, 4t) = / p(y, 4t) rfy 



Trr^ Jr 



(and wo{r) :— 8r ^uo(f^)), problem (|2.ip - p.4p becomes equivalent to 



wt — Aw 



';WWr, < r < 1, t > 0, 



Wr{0,t) =^ 0, i>0, 
w(l,i) = 8, 00, 
w(r, 0) = wo(^)i < r < 1, 



(3.17) 
(3.18) 
(3.19) 
(3.20) 

:Wr, which in turn corresponds to the radial Laplacian in 4 space dimensions. It 
should be noticed that w has the same scale invariance as p, but is smoother than p. Problem p.l7p - 
(|3.20p turns out to be convenient regarding the study of C^ regularity of u up to the boundary, which 
was left open in [5] . Namely, using this transformation, one can show the following additional properties 
for u. 



where Aw 



Wr 



Lemma 3.2. Let uq satisfy 
For any given T > 0, if 



then 
and 



and 



and denote by u the global solution of problem I12.1]) - (K^ . 
u(x, t) 



sup 



< 00, 



0<a;<l, 0<t<T X 

u, eC([0,l]x(0,T]) 

u^(o,i)>o, o<t<r. 



(3.21) 

(3.22) 
(3.23) 
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Note that the assumption (|3.2ip (for all T > 0) will be shown in Section 4 (see Lemma 4.3(ii)), as a 
consequence of our main upper solution construction, hence leading to the global C^ regularity property 
(Pl^ in Theorem O 

Proof of Lemma \S.2[ Since the semilinear parabolic equation in (|3.17[) has only linear growth with 
respect to the gradient, standard arguments based on the variation-of-constants formula (see, e.g., [431 
Example 51.30], or [151 P- 889]) show that problem (|3.17p - (|3.20p is locally well-posed in the space of 
(radial, nonncgativc) L°° functions. More precisely, for any < wq G L°°{0,1), there exists a unique, 
maximal (radial, nonncgativc) classical solution w of p.l7p - (|3.20p . with w S C^'^([0, 1] x (0,Tm)) and 
w e C([0,T™);L9(0, 1)) for aU finite q < 1. Moreover, 

Tm < oo => lim ||iu(t)||oo ~ oo. (3.24) 

Now, if uq satisfies (|2.5p - (|2.6p . then woir) := 8r^^UQ{r^) verifies < wq £ L°°{0, 1). Denote by w the 
corresponding maximal solution of (|3.17p - (|3.20p and let 

u(x, t) = ^w^yf^, t/^), 0<a;<l, 0<t< AT^. 

8 

We see that u satisfies the regularity conditions in (|3.8p . (|3.9p and (|3.1ip with r = 4T,„. Also, since 

1 fx 

u^(a;,i) = ^w(^/^,f/4)^-^u.,(^/i,t/4), < x < 1, < i < 4T„, 
o lb 

we have 

weCi'"([0,l] X (0,4T„)), (3.25) 

hence in particular (|3.10p with r = 4T„i. Then one easily checks that u solves p.ip - ()2.3p on [0, 4Tm), along 
with (j3.7p . We may thus apply Proposition 13. II to deduce that m = u on (0,To) with To = min(r, 4Tm). 
We claim that T™ > T/4. Indeed, if Tm < T/4, then ([3:241) imphes 

u{x,t) 
lim sup = oo, 

contradicting ([3.2ip . Consequently, property ([3.22p follows from ([3.25p . Moreover, w > in [0, 1] x (0, T„i) 
by the strong maximum principle, which readily implies (|3.23p . D 

3.3. Small time estimates. Throughout the paper, we denote by V the parabolic operator defined by 

Vv := vt - xvxx - '^vvx- (3.26) 

The following two lemmas will be useful to initialize the comparison between u and the main lower/upper 
solutions constructed in Section 4. 

Lemma 3.3. Let uq satisfy 112. 5]) and V2. 6\) and denote by u the global solution of problem, V2.1\) - ^K^ . 
Then there exist t^tj > such that 

u{x, t) < 2Kx, < a; < 1, < t < r, (3.27) 



and 



Proof. Define 



Since 



u(a;,r) < l-?7(l-a;), < x < 1. 



v{x,t)= ^^^^^ , 0<x<l,0<t<l/2K. 

2K^x 2K^x 

Vv = = 

(1 - 2A'i)2 (1 - 2Kty 



and u(l,t) > K > 1, the comparison principle guarantees that 

u{x, t) < v{x, t) < 2Kx, < a; < 1, < f < 1/AK. 

Fix T = l/4_ftr. By Hopf's lemma, we have Ux{\,t) > 0. Since m(-,t) is nondecreasing in x, this implies 
u{x^ r) < 1 — 7/(1 — a:), < x < 1, for 77 > sufficiently small and the conclusion follows. n 
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Lemma 3.4. Let uq satisfy i2. 5\) and h2. 6\) and denote by u the global solution of problem i2.1]) - (K4^ . 
For any given 6 G (0, 1), there exists Ts > such that 

u{x,Ts)>mm{l-6,x/5), < a; < 1. 

Proof. By Lemma 13.21 we may fix a small r > such that Ux{-,t) e C([0,1]) and Ux{0,t) > 0. 
Therefore u{x,t) > ijx for all x G [0,1] and some ?/ e (0,1)- Since V[rix] < and ?; < 1 = u{l,t), the 
comparison principle implies that 

u[x,t)>'qx, Q<x<l,t>T. (3.28) 

On the other hand, by (|2.7p . there exists T > t such that 

u{x, t) > 1 - (5, il-S)S <x<l, t>T. (3.29) 

Define 

v{x,t) = {7] + 2if{t-T))x, 0<x<l,t>T. (3.30) 

We have 

Vv = 2rj'^x - 2(?7 + 2jf{t - T)fx < (3.31) 

and, due to ((3?28)) . 

u{x,T)>v{x,T), 0<x<l. (3.32) 

Since rj < 1 < 1/6, we may find Ts > T such that 

i1 + 2tj\Ts-T) = 1/5. (3.33) 

In view of (0221), (^3(1)) and ([3^, we have 

u{{l~5)5,t) >l-5>v{{l-S)5,t), T<t<Ts. (3.34) 

It then follows from p.3ip . (|3.32p . p.34p and the comparison principle that 

u{x, Ts) > x/S, < a; < (1 - 6)6. 
Using 'P[x/6] < and ([X^ . we deduce that 

u{x, t) > x/6, < a: < (1 - 6)6, t > Ts. 
This combined with (|3.29p yields the desired conclusion. D 

3.4. ODE lemmas. We consider the differential operator 

„ 2yw' 2w 

Lw := yw +- \- — rrr. 

1 + y (1 + 2/)2 
First, the expression for the operator C reads 

^y + lJ V y J 
If ip E C{[0, oo)) satisfies ipiy) = 0{y) as y ^ 0, then the problem 

Cw ~ %Ij, y > 0, 
w(0) =0, w'{0) = 
admits a unique solution w, and w can be represented as 

(note that the integral is convergent due to the assumption on ^jj). In particular, 

V- > on [0, oo) =^ £0 V' > on [0, 00). (3.35) 



Cw 
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On the other hand, wo{y) = / ^lyi solves 

Cw^O, y > 0, 

w{0) ==0, w'(0) = 1. 

In the next section, to construct our main upper and lower solutions, we will need to know the 
asymptotic behaviour of the action of the operator Co on some particular functions as y -^ oo. More 
precisely, let us define 

^-uf y 



and 

where M > and 



/ = 2//' - yf + /, 
9 = ^o'f 



h = £-!(/ + M^) =g + MC^'ifi, 



^eCl([0,Oo)), ^(0)^0, ^(y)= y>2 



logy 



We have the following lemmas. 

Lemma 3.5. As y — > cx), the function f satisfies 



f{y)=logy-2 + 



(i) 
(ii) 

y "• y 

Lemma 3.6. As y -^ oo, the functions g and h satisfy 

y^ogy _ 9y 
2 4 



log^y 



(i) 
(ii) 
(in) 



liy) 



0{\og^y) 



W2/) = ^--4+0(^ 



Kv) 



ylogy 9y 



4 ^""UogyJ' 



(iv) 



^"' 2 4 Vlogy 



Proof of Lemma \3.5\ (i) Usin 



log(y + l) ==logt/ + 0(l/y), asy^oo, 



we obtain 



/(y) 



(j/ + l)2 

y 


(1 + 2/)^ 

y 



t 



2 „t 



Vs + 1 {s + lf 



ds dt 



^ ^1 + ^ + ^) (log(t+l) 



t + 1 



dt 



[i + y? 



[ylogy - 2y + 0{\og^ y) 



logy -2 + O 



log y 



as 2/ — > (X). 



(3.36) 
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(ii) Now using (i) we deduce 

f'{y) - \--^ 
[y y + i 



/(y) + i(^log(, + l)-^ 



I f\og^y\ 

= — h C — 5— , as 2/ ^ 00. 

y \ y J 



To show Lemma 13.61 we first note that, due to Lemma 13.51 and (j3.36[) . we have 

f{y) = /i(y) - 3/2(2/) + 0{h{y)), as 2/ ^ cx), 
where 



/i(2/)=log(2/ + l), /3(2/) 



y + 1 



and 



/2eCi([0,oo)), /2(0) = 0, /2(y) = 1, 2/ > 1. 

Denote pi = £q /i, for i = 1, 2, 3, and 54 = £q (p. Lemma 13.61 wih then be an immediate consequence of: 
Lemma 3.7. As y -^ 00, the functions gi satisfy 



(i) 
(ii) 

(iii) 
(iv) 



ylogy Sy , . , logy 1 ^f^ogy 
5i(y) = ^ -j + 0{logy) and g-^{y) = — ^ + 0(^ 



92{y)^l + 0{l) and g'^{y)^\ + o(- 
2 2 \y 



gM^0(\0g^y) and g^(y) = o(l^^J^) ; 



Proof, (i) As 2/ — *■ 00, we have 
51 (y) 



(2/+ 1)^70 V t 



y ft+\ 



j ((t + l)log(f+l)-t)di 



{i + yYJo 



[(i + 1) log(t + l)-t + {\og{t +l))\dt 



y 



y 



2/2 log y 3j/2 



0(2/ log 2/) 



hence 



9'M 



1 



2/ y + 1 

-i + 

2/ 



2/ log 2/ % , ^n ^ 



1 /"^ 
gi{y) + - / log(s + l)ds 



2/ log 2/ 32/ 



+ 0(log2;) + 



(2/+l)log(2/ + l)-2/ 



log 2/ _ i_^o / log 2/ 
2 4 \ y 
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(ii) Due to the definition of g2 there exist constants Ci, C2 G M such that, for y > 1, 

2 
(Ci + t)dt 



52 (y) = 

Therefore, for y > 1, 

92{y) 



(iii) As y — > cx), we have 

53(2/) 



y 



(2/ + 1)2 

y 

(2/ +1)2 [ 



^ n + iV r* 

f2{s)dsdt + 

^1 



t 



t 



1 


2 


2/ 
1 


2/ + 1. 
2 


2/ 


2/ + 1. 



C2 + 1- + 0(2/) 



.92(2/) 



^ + 0(1), asy^oo. 



(2^ + 1)^ 



2/ + 1 

y 



[Ci + y) 



.g2(2/) + l + 0( ij =i + o(^i ) , asy^oo. 



log'''(t + l)di 



3(2/ +1)^70 V i 



y 



0{\og^{t + l))dt 



m+yrJo 

^^ly^2 O ((2/ + 1) log' (2/ + 1)) = O {\oi y) , 



hence 



53(2/) 



1 2 

2/ 2/ + 1 
1 



1 [y iog^^i+2) ,, 

53(2/) + - / -^ -^dt 



, 0[-,))Oi}oey), 

(iv) Similarly to (ii) in this case we have, for y > 2, 



1 



yJa (i + 1) 

\og\y + \) 



32/ 



O 



log y 



9i{y) 



y 



(2/ + l)^ 


y 


(2/ + 1)2 


2/ 



2 /i+lX2 .* 







t 



ip{s)dsdt - 
J2 



t+1 



Ci 



ds 
logs 



dt 



-':(T)°fe 



di 



ol'^l^o 



(2/ + 1)2 Vlog2/ 
Therefore, for y > 2, 

^1 2 



2/ 



logy 



as y — *■ 00. 



94iy) 



y 2/+ 1 



.94 (y) 



y 



(y + i)^ 



y + i 



Ci 



y ds 



logs 



o 



logy 



as y -^ cx). 



D 



4. Proof of the main results 

4.1. Construction of a lower solution. Motivated by the idea of an asymptotic expansion around 
(moving) steady states and based on a self-similar variable, see Subsection 1.3, we make the following 
Ansatz: 

u{x, t) = l- -^ + b[t)f{y) - b^{t)g{y), y = a{t)x, 

y + 1 

where the functions a, b, /, g have to be determined. Note that the variable y now ranges into the time- 
dependent interval [0, a(t)]. Here, a and b are expected to satisfy 

a{t) ^ Ux{0,t), hm a{t) = 00 and hm b{t) = 0. 
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Lemma 4.1. The problem 



Ik ~ ^Ikxx ^ SmUj., 0<a::<l, f>0, 
u(0,i) = 0, t>0, 
u(l,t) < 1, t>0, 



admits a solution of the form 

u{x, t) = 1 - -^ + h{t)f{y) - b^{t)g{y), y - a{t)x, 

y + I 

where the smooth functions a>0,b>0, f>0 and g have the following properties: 

f{y}^\ogy, g{y) — , as y ^ oo, 



(4.1) 



/(0)=5(0) = /'(0)=5'(0) = 0, 
"5 



(t) = (l + 0(t-^/^\ogt))cxp 



'2t 



, as t ~y oo, 



ana moreover 



^^^ l + 0(f-l/2) 

bit) = — -— ——, as t — > (X), 

^ ^ a{t) log a{t) 



u^>0, < .T < 1, t > 0. 



(4.2) 
(4.3) 

(4.4) 
(4.5) 



Proof. Step 1: Construction of the subsolution. In what follows we shall omit the variables t 
and/or y when no confusion is likely. We take u as in (|4.ip where we assume 

(4.6) 
(4.7) 



a > and lim a{t) = oo. 

t — ^oo 



We compute 



ii + y? 

-2a2 



baf'{y)-h'ag'{y), 



and 



hence 



a' y 



,,^ ,,+ba'ny)^b'a'g"iy), 
(1 + yr 

b'fiy) + ba'xfiy) - 2bb'g{y) - b'a'xg'{y) 



b'fiv) + —yf'iy) 2bb'g{y) - —yg'{y) 
a a 



^'" a{l + yr 
Recall that the operator V is defined in ()3.26p . It follows that 



'Pu = 



a y 

a (l + 2/)2 

+ 2a- 



y 



b'f + —yf - 2bb'g - — yg' 
a ' ' a ' ' . 



y 



{i+yr 

Collecting terms of same order in b yields 



- ba yf" + b^a yg" - 2a -^^— + bf - Vg 



i + y 



1 



+ fc/'-6V 



Vu 



_a y 

a [l + yf 



J -\ yj -2bbg yg 

a a 



ab 



yf" 



{l + yr 

, 2/ , 2yf' 



(l + y)2 l + y 



-ab" 



yg 



2g , 2yg' 



{1 + y)^ l + y 
The natural scaling of the equation leads to the choice 

a' 



2ff' 



2ab- 



f'g + fg' 



4„„' 



2ab gg 



(4.8) 



(4.9) 
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SO that in the RHS of 
Assuming 

we also denote 
and observe that 



where the last equality conies from 



the first term will be of the same order as the terms in the second bracket. 



a' > 0, 



7:= - 



- = ab, b' = -{l + j)ab^, 



(4.10) 
(4.11) 

(4.12) 



a\' / 1 \' 



T=U7 =[7rh 



1. 



■.a'/ \abJ ab^ a^b ab^ 

To make things clear, let us already stress that the final choice of a will guarantee 

7 > and lim j{t) = 0. 



Using (14.121) . identity (|4.8p can be recast under the form: 

y ...„ 2/ 2yf 



Vu 



ab 



[i + yy 



vl" 



{l + vf l + y 



+ab^ 
+ab^ 



yg 



'' +^-2//' + ,/' -(1+7)/ 



(i + yy- 'l + y 

2/'.g + 2/.g'-yg' + 2(l+7)5 



4„„' 



2ab-'gg^ 



hence 



Let us now choose 



Vu = ab 



(l + y)2 
+ab^ 



Cf 



ab-" 



Cg-2f,r + yf'-{l+j)f 



2/'.g + 2/5' - yg' + 2(1 + 7)5 - 2ab^gg 



f:={I + C^') 



y 



> 



y 



>o, 



(4.13) 



(4.14) 



,(i + y)V - i^ + y? 

which solves Cf = y/(l + y)^ for y > 0, with /(O) = and /'(O) = 1 (cf. Subsection 3.4). Dividing by 
ab^ , we see that Vu has the same sign as the quantity 



Next we choose 
Therefore 



Cg - 2//' + yf - (1 + 7)/] + b [2fg + 2fg' - yg' + 2(1 + ^)g\ - 2b^gg'. 
g:=C^\2fJ'-yf' + j). 



A = -jf + b 



2fg + 2fg'-yg' + 2{l+j)g 



2„„' 



2b'gg 



We now proceed to show that the quantity A is nonpositivc for large t, by considering separately the 
regions yo < y < a{t) and Q < y < yo, for some large j/o independent of t. At this point, we make the 
additional assumptions that 

1 

7 ^ ., _ ^ , as i — *■ 00 (4-15) 



and 



log a' 



hence b ■ 



1 



-, as t — > 00 



, , , , „. . ^„ (4.16) 

log a a log a 

(which will be verified on the final choice of the function a - actually more precise expansions of 7 and b 

will be needed in the final matching process) . By Lemmas 13.51 and 13. 6i we have 

/~logy, f'g^o{g), Jg' = o{g) and yg' -^ g — , as y ^ 00. 
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Consequently, fixing S > and taking j/o and to large enough, we have, for yo < y < a(t) and t > to, 

-7/<(-l + <5)|^^, g,g'>0, 
log a 

2f'g + 2fg' - yg' + 2(1 + 7)3 <5y\ogy- -y\ogy + 2(1 + 7)(2 + 5)y\ogy < (3 + ^^)y^ogy. 



and 



Taking 6 — 1/12 and also using (|4.16p . we thus obtain 

A<[-l + 6)'^ + ^(\ + ,6)ylogyJ^(-l^5+(\ + ,6y-)<, 
log a a log a V 2 / log a\ V 2 / a / 

for yo < y < a(i) and t > to (possibly larger). Next, for < y < yo, (|4.14p implies 

/(y) > ciy, Q<y <yo, 

whereas /(O) = g{0) = yields 

|2/'g + 2fg' - yg'\ + 4|.g| + \gg'\ < c^y, 0<y<yo, 

with ci,C2 > 0. Consequently, 

A < -ci7(t)y + C2bit)y = [-ci + ca-^l 7(^)2/ < 

L 7(t)J 

on [0, yo] for i large enough, due to (|4.15p and (|4.16p . 

We have thus proved that, under conditions (|4.6p . (|4.10p . (|4.15p and (|4.16p . there holds Vu < in 
(0, 1) X (T, 00) for T large enough. By a time-shift we may obviously take T — 0. 

Step 2. Determination of a{t) by matching at the outer boundary. Next, the determination 
of a(t), will be done by "matching" with the boundary condition at x = 1, i.e. by writing 

u{l,t)<l, 

which is equivalent to 

bfia)~b^gia)<^. (4.17) 

a + 1 

It is of course sufficient to check (|4.17p for large t (thanks to the possibility of shifting time). Let us 

first sketch the resolution of (|4.17p in a rough way. Using (|4.9p and applying Lemmas IS.Sf i) and I3.6f i) at 

leading order, we are left with 

a' /, o-' a log a \ ^ 1 ,, ,„s 

^(loga-^^j<^, ast^oo. (4.18) 

We expect the second term in the bracket of the LHS of the preceding relation to be much smaller than 
the first one as a —> 00. If we ignore it we obtain the differential inequality 

a'<-^, asi^oo, (4.19) 

log a 

which implies that 

ait) < e^\ as t -^ 00. 

However, the latter estimation is not accurate enough to show the desired estimate and we thus add a 
correction term in (|4.19p . More precisely, we look for a{t) as the solution of 

a' = -^(1 + 77), t>0, a(0) = 2, (4.20) 

log a 

where the correction term 77 ~ 77(a) has the form 

7/=—^ + —^, /v>0. (4.21) 

21oga log^a ^ ' 
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Now plugging (|4.20p into (|4.17p . recalling also (|4.9p . and using the exact asymptotic behaviour of 
f{a),g(a) a,s a -^ oo, provided by Lemmas 13. Sf i) and 13.61^ 1). we are reduced to the condition 



1+7? 

a log a 



^ 1 log a\ 1 + ri f a log a 9a ^ /, ■? 



< 



1 



a + l' 



as t ^ oo 



or 



(1 + ^/) 



5 + 7] 



9r/ 



2 log a 4 log^ a 4 log^ a 
Plugging ((42T|) into ((422)) we first obtain 



O 



K 



2 log a log^ a 



1 



and finally 



^ 4is: - 21 ^ 

1+ 2— + 

4 log a 



2 log a log a 
1 



O 



log a 



1 



< 



a + l 



as t -^- oo. 



(4.22) 



log^a 



< 1- 



log'^a 
1 



< 1 



a + l' 



as t ^ oo 



a+V 



as i — > oo. 



(4.23) 



In order for (|4.23p to be satisfied, we choose K < ^ and we then obtain the following ODE for a: 



log a 



fl + — ^ + — ^y OO, a(0) = 2. 

V 21oga loe^ay ^ ' 



(4.24) 



Here we should mention that the same form for 77, given by (|4.2ip . will be also considered for the upper 
solution, with a different constant K; see next subsection. Equation (|4.24p implies that 

logo -5/2 

and integrating with respect to t we get 



(l + 0(log a)) , as t ^ 00, 



log^ a - 5 log a = 2t + O 



ds 



as t ^ 00. 



/o log {a{s)), 
Solving the quadratic polynomial in log a, and noting that log a(s) > \/2s by (|4.24p . we end up with 



Tm. As for ga, it follows from 



fo= ^ = 



K \ 1 + (log-i a) _ 1 + O (i-i/2) 



a? a log a \ 2 log a log^ a 



a log a 



a log a 



as i ^ 00, 



(4.25) 



where we used (|4.24[) and log a > v2t, as t -^ 00. On the other hand, denoting G{s) = s + (5/2)s^ + ifs 
and using (|4.24p . wc sec that 7 = {a/a'Y satisfies 



1 



G'(l/loga) 



G'(l/loga) 



= H(l/loga), 



(4.26) 



G(l/loga)J alog^a G2(l/loga) log^ aG(l/loga) 

where H{s) = s(l + 5s + 3^s2)(l + (5/2)s + A's^)"^ Finally, the assumed properties (|46)) . (|4+0)) . (|4+5l) 
and (|4.16p of a, b, 7 are immediate consequences of (|4.24p . ()4.25p and (|4.26p . 



Step 3. Proof of (jiTSJ) . By (gj]) we have 



1 



a u^ = 



+ 6/'(y)-5V(y). 



-" {i + vY 

By Lemma 13.61 and ()4.16p . taking yi and ti large enough, we have, for t >ti and yi < y < ci{t), 



a ^u^ > 



1 logy ^J_/l 



1 



2a2 a2 log^ a a? \2 log a 
Now, for < J/ < 2/1 and t>ti possibly larger, we get 

1 



>0. 



a u^> 



Cb{t)^Ch\t) >0. 



By a time-shift we may obviously take fi = and (|4.5p is proved. 
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Remark 4.1. It is still possible to obtain a qualitatively correct lower solution with just a two-term 
expansion, for instance by making the simple choice f{y) = log{l + y), <? = in (|4.ip . However this yields 
only the lower grow-up rate up to a multiplicative constant, i.e. Ux{0,t) > Ce^^*, and docs not enable 
one to deduce an expansion of the form 



4.2. Construction of an upper solution. The form (14. ip does not seem to be sufficient to construct 
an accurate upper solution (i.e., leading to a function a{t) fulfilling (j4.3p ). We need a slight perturbation, 
corresponding to the modified Ansatz: 

u{x,t) = l-^+b{t)f{y)-b\t)g{y,t), y = a{t)x, (4.27) 



where 

and e{t) goes to as t ^ oo. 
Lemma 4.2. The problem 



y+l 

~g{y,t) = {l+e{t))hiy), 



Ut — XUxx > 2MM3;, < X < I, t > 0, 

u{0,t) ^ 0, <>0, 
u(l,t) > 1, t>0. 



(4.28) 



admits a solution of the form, 
u{x, t) = 1 



y+l 



+ b{t)f{y) - (1 + e{t))b\t)h{y), y = a{t)x, 



where the smooth functions a{t), b{t), /(y), hly) have the following properties: 

2/ logy 



f{y)^\ogy, h{y) 



as y —f 00, 



/(O) - h{Q) = /'(O) = h'{Q) = 0, 
5 



i{t) = h +0(f"i/2logt)jexp 



as t ^f 00, 



as i —> 00, 



(4.29) 



(4.30) 
(4.31) 



a(t)loga(i) 
e(t) ~ (2i)"^/^ as t^ 00. 

Proof. Step 1: Construction of the supersolution. Taking u as defined in (|4.27p . (|4.28p . the 
expression for Vu is similar to (|4.13p . except that g, g' , g" arc now replaced with 5, gy,gyy and an additional 
term —b^e'h is added, which is inherited from ut- As in the proof of Lemma 14.11 b and 7 are defined 
through (|49l) and (gjH), and we assume ([46ll . (|4?T0l) . ([4T5| and ((4T6ll . This leads to 



Vu 



ab 



{1 + yy 



Cf 



ab" 



Cg-2ff + yr-{l+j)f 



-a¥ 



2f'g + 2f~gy - ygy + 2(1 + 7)3 



-2ab^- 



Replacing g{y,t) with (1 + e{t))h{y), wc obtain 



Vu^ 



ab 



{i + yy 



Cf 



ab-" 



{l + e)Ch-2ff + yf'-{l + j)f 



Denote 



+ {l + e)ah^ 2f'h + 2fh' - yh' + 2(1 + -f)h - 2(1 + efaf^hh'. 
Bo := (1 + e)Ch - 2//' + yf - (1 + 7)/ - -h. 



99y 



(4.32) 
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As in Lemma HTTl we first choose 

Dividing identity (j4.32p by ab^, we see that Vu has the same sign as the quantity 

B := Ba + {1 + e)b\2f'h + 2fh' - yh' + 2(1 + -/)h\ - 2(1 + eflPhh'. 

Next we choose 

h^Q\2fJ'-yf' + f + Mv), 
where Lp satisfies 

ifiQ) == 0, ip'{0) > 0, Lp{y) > for < y < 2, Lp{y) = 1/ logy for y>2. 

Note that, taking M > 2 suitably large, we have 

2ff'~yf' + f + Mip>0, y>0, 

by Lemma |3.5[ hence 

Hy) > 0, 2/ > 0, 
due to (|3.35|) . We compute 

So = (l + £)(2//'-y/' + / + M^)-2//' + y/'-(l + 7)/--/i 

a 

= e(2//' - yf) + M(l + e)^ + (e - 7)/ - -h. 



(4.33) 



(4.34) 



At this point, we choose 



£ = 7, 



(4.35) 
(4.36) 



where 7 is defined by (|4.1ip . and we assume again (|4.15p and (|4.16p . along with 

7' <0 
(these assumptions will be verified on the final choice of the function a). By (|4.34p and (|4.36p we have 

Bo>l{2ff'-yf')+Mil+j)ip 
and it follows that 

(1 + j)-^B > \—^{2ff' - yf) + mJ +b\2f'h + 2fh' - yh' + 2(1 + j)h\ - 2(1 + j)b^hh'. (4.37) 

LI + 7 ' J L J 

To show that the RHS of (|4.37p is nonnegative for large t, we again consider separately the regions 
yo l£ y l£ o-it) a-iid < y < yo, for some large yo independent of t. By the estimates in Lemma 13. 5[ we 
have 

f^logy, f'^- and yh'^h^^^^, as y ^ 00. (4.38) 

y 2 

Consequently, fixing S > 0, using (|4.15p and taking yg and tg large enough we have, for yo < y < a{t) 
and t > to. 



-^(2//'-y/') + M^>-|7+T^> ^ 
1+7 2 logy 



M M - 2 
> 



3gy log a logy logy 

2f'h + 2fh' - yh' + 2{l+j)h >-{l + s)ylogy + 2(1 + 7)(^ - s)ylogy > (i - 4s)ylog, 



and 



2(l + 7)M'<ylog2y. 
Assuming M > 3, taking 5 = 1/8, and also using (|4.16p . we infer that 

N 1 „ 1 ylog y 1 / ylog y\ If loga^, 



logo a^log^a logo V a^loga/ logo 
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for Ho < y < a{t) and t > to. Next, for < y < yoi (|4.33p implies 

M(p{y) > ay, 
whereas /(O) = h(Q) = yields 

2ff'-yf'>-C2y, 2f'h + 2fh'-yh' + 2{l + j)h>~C2y, 2(1 + j)hh' < C2y, 
with ci,C2 > 0. Therefore, 



(1 + j)-^B > -jit)c2y + ciy - h{t)c2y - b'{t)c2y = ci - 02(7^ + h{t) + h'it)) 



2/>0 



on [0, j/o] for i large enough. 

We have thus proved that Vu > in (0, 1) x (T, oo) for T large enough. By a time-shift we may 
obviously take T = 0. 

Step 2. Determination of a(i) by matching at the outer boundary. The determination of a{t), 
will be done again by "matching" with the boundary condition at a; = 1. Imposing that 

u(l,t) > 1 

we obtain 

1 



hj{a)-¥{l + e)h{a)> 



a + 1' 



and taking (|4.9p into account we end up with 



-[f{a)--{l + e)h{a)] > ^. 



(4.39) 



Again it suffices to check (|4.39p for large t. Following the same reasoning as in the case of lower solution 
we again look for a{t) as a solution of 

a 



log a 



(1 + ?7), as < ^ 00 



(4.40) 



where the correction term 77 = 77(a) is given by (|4.2ip with K a constant to be determined. 

Plugging (|4.40p into (|4.39p and using the exact asymptotic behaviour of f{a),h{a) as a — > 00, given 
by Lemmas 13.51 and 13. 6i we arrive at 

'log^a\ (l + 77)(l + e) /a log a 9a 



1 + 7? 

a log a 



loga-2 + O 



or 



(1 + ^) 



1 



O 



a log a 
loga\ (l + ?7)(l+e) /I 



2 4 ^ Vloga 



> 



a+1 



9 



O 



f 1 



> 



a + 1 



(4.41) 



log a ' \ a J log a ^ 2 4 log a ' \ log^ a , 

as i ^ 00. Denote by T the quantity in the bracket in the LHS of (|4.4ip . Using (|4.2ip and (|4.4ip we 
obtain 



^_ 2 ^^^loga 
log a 



K 



2\oga ' log^aj \' ' loga ' 21og^ay V^foga 41og^a " Vfog^a 



' o( ' 



o 



f 1 



2 log a 2 log a Vfog 

Then (|4.4ip becomes equivalent to 
5 K 



as i — > 00. 



2 log a log-^ 



2 log a 2 log 



O I 5 — I I > 1 , as i -^ 00 

Vlog^ 



that is. 



4/\ - 23 ^ / 1 \ 1 

1 H T h O — o — > 1 7, as i ^ 00. 

41og2a Vfog a/ ~ a+1 



(4.42) 
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For (|4.42p to be satisfied we choose K > ^, and we again take a to be the sohition of the ODE 
(g^. Then, by the end of Step 2 of the proof of Lemma [ill we obtain gJH), ((430)l . as well as the 
assumed properties gj]), ([iTTO]) . (|4?T5l) and ([iTTe]) of a, 6,7. Finally, we note that (|4?T5l) . ([429]) and 
(|435l) guarantee ([43T|) . and that ([426]) imphes 

7' = — ^i/'(l/loga) = — i-(F'G)(l/loga) ^ — ^, as i ^ 00, 
a log a log a log a 

hence (|4.36p (after a further time-shift). D 

4.3. Proofs of Theorem [27T] and Corollary [2T2I Let u be the solution of ((2?T|) - ((2^ and let u,u he 

the lower/upper solutions provided by Lemmas 14.11 and 14.21 The asymptotic expansion (|2.8|) - ()2.9p in 
Theorem 12.11 will be an immediate consequence of the following two Lemmas. The first one guarantees 
that u lies between suitable time-shifts of u and u. The second one shows that (shifted versions of) u 
and u satisfy the required asymptotic behaviour. 

Lemma 4.3. (i) There exists Ti > such that 

ui-,t)>u{-,t-Ti), t>Ti. (4.43) 

(ii) Let T be as in Lemma \3.3[ Then there exists T2 > such that 

u{-,t)<u{-,t + T2), t>T. 

Proof, (i) Since u(-,0) € C\[0,1]) with u(0,0) = 0, u(l,0) < 1 and u^(-,0) > 0, it follows from 
Lemma [3.41 that u{-,Ti) > w(-,0) for some Ti > 0. The assertion then follows from the comparison 
principle. 

(ii) Due to (ji^S]) . (|4T4l) and h{0) = 0, we have 

/(y) > ci log(j/ + 1) and /i(j/) < C2 (y + 1) log(y + 1), y > 0, 

for some ci, C2 > 0. This along with (|4.29p - (|4.3ip implies that, for all t > t2 large enough and all x G [0, 1], 

u{x, t) = —^ h bf{ax) - (1 + e{t))h^h{ax) 

ax + 1 



ax 



ax 



- + &log(aa; + l)(ci - 2c2(aa; + 1)6) 



ax , , , , s / „ a + 1 \ CLX 

> hbloe ax + 1) ci -3c,— > 



— + 61og(aa; + l)(ci - 3c2— ^ ) 

1 \ a log a / 



ax + 1 V a log a J ax + I 

Take t, rj as in Lemma [373l and set xq = l/AK. For all x e [0, xq] and t > t2, with t2 > t possibly larger, 
we have a/2K > ax + 1, hence 

QjX 

u(x, t) > > 2Kx > u(x, r), < x < xq, t > ^2- 

ax + 1 

On the other hand, (|4.38p implies that 

./'(y)<C3 and |/i'(2/)|<C3 + log(2/ + l), y>0, 

for some C3 > 0. For all x £ {xq, 1] and t > t2 (possibly larger), we thus have 

+ baf'{ax) - (1 + £{t))b^ah'{ax) < — — + c^ba + 2{cz + log(a + l))b^a < 77, 



(l + aa;)2 •' ^ ^ v wy v y- (1 + a2.Q)2 

due to ([i?^ and (|T^ . hence 

u{x, t) > 1 — 77(1 — x) > u{x, r), Xq < X < I, t > t2. 

Therefore u{-,t2) > u(-,r) in [0, 1]. The assertion, with T2 = ^2 ~ ■'', thus follows from the comparison 
principle. n 
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Lemma 4.4. Let 



Then, for any T G M, each of the functions w ~ u and w = u satisfies 

1 



A{t) = cxp 
ons 
l~w{x,t + T) = 



1 + A{t)x 



l-x + 0(t~^^'^\ogt) 



as t ^ oo, uniformly in [0, 1]. 



Proof. Wc shall give the proof for w ^ u, the other case being completely similar. Set d{t) = a{t + T) 
and b{t) = b{t + T). We first note that, by ([43ll . we have 



~a{t)={l + 0{t-^'^\ogt)^A{t) 



and 



\og~a{t) = log A{t) + 0{t-^l^ logt) = (^1 + 0{t-^ logt) j log A(t). 
Also, by (|4.4[) . we have 

a log a 
Moreover, due to Lemma [3.51 there exists C > such that 

|/(y)-log(l+2;)|<C, y > 0. 

Now, using (|4.45p . (|4.46p and (|4.44p . we compute 

1 r., . 1 



(4.44) 
(4.45) 
(4.46) 



1 + aa; 



hf(ax) 



1 + ax 

1 
1 + aa; 

1 
1 + aa; 

1 
1 + a.a; 



^ _ (1 + ~ax)f{hx) ^ ^ ^ (5^^-1/2 



1- 
1- 



aloga 
(1 + aa;)log(l + aa;) 

a log a 
a;log(l + aa;) 



l + 0(i-l/2^ 



where 



R{x,t) 



, . +0(t-i/2) 

log a 

l-x + R{x,t) + 0{t-^/'^) 

a; (log a — log(l + da;)) 
logd 



Here and in what follows, the O's are uniform in [0, 1]. To control R we note that, if da; > 1, then 

log(da;) < log(l + da;) < log(da;) + log 2, 

hence 

.(|log.|+log2)^ C 

log a log a 

whereas, if ax < 1, then |i?(a;,t)| < a; < 1/d. It follows that sup^^grQ -n |i?(a;,t)| = 0(f^^/^); as i — > oo. 
Since, by (|4.3|) . we have 

0(t-i/2 logt)^a;' 



1 + da- = (1 + ^a;) 
we deduce that 

hence 



(d - A)x 



1 + Ax 



(1 + ^a;) 



1 



1 



1 + da; 



6/ (da;) = 



1 + Ax 
l + 0(t~^l'^\ogt 



(l + Aa;)(l + 0(t-i/2logf)), 



(l_^ + 0(t-i/2)j, 



1 

1 + da; 



If {ax) 



\ + Ax 

_ l-a; + 0(t-i/2log<) 
~ 1 + ^a; 



(4.47) 
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On the other hand, due to Lemma \W^ and g{0) ~ 0, there exists Ci > sueh that 

l5(2/)|<Ci(l + 2/)log(l + 2/), y>0. 
Consequently, using also (|4.45p . we obtain, for t ^ oo, 

(1 + A.)b'\gi~a.)\ < 2Ci (^ + «^)^ l"f + «^) < ^ ^ o(,-i/2). 

a^ log a log a 

Combining (|4.47p and (|4.48p finally yields 

1 . ^ ^ -.. l-x + Oft^i/^logt) 

l-u{x,t + T) = - —-bf{hx)+b^g{dx) ^ ^' 



(4.48) 



1 + ax 1 + Ax 

□ 

Proof of 112.1 (J\) and V2.ll]) . First notice that estimate p.27p in Lemma [5751 and Lemma HT^T ii) guarantee 
the control of the slope at a: = 1, namely (|3.22p . for any finite T > 0. The C^ regularity property (|2.10p 
is then a consequence of Lemma F3. 2 1 To show p. lip , note that M^(0,i) = a{t) due to (|4.7p and (|4.2p . 
The lower estimate corresponding to ()2.1ip is then a consequence of (|4.43p and ()4.3p . The proof of the 
upper part is similar by using u. D 

Proof of Corollary [KM Assertion (i) is an immediate consequence of p.Sp and (j2.9p . To show (ii), it 
suffices to observe that, due to (|2.9p . 

1 



dx 



/o 1 + A{t)x 
as t ^ oo, and to use 



\og{l + Ait)x) 
Ait) 





1 



Mi±p.(i+or./=))V2< 



exp 



-l-VTt 



X dx 1 

< 



1 + A{t)x - A{t)' 
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